Wigner flow reveals non-classical features of quantum phase space dynamics 
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The behaviour of classical mechanical systems is characterised by their phase portraits, the col- 
lections of their trajectories. Heisenberg's uncertainty principle precludes the existence of sharply 
defined trajectories, which is why traditionally only the time evolution of wave functions is studied 
in quantum dynamics. These studies are quite insensitive to the underlying structure of quantum 
phase space dynamics. We identify the flow that is the quantum analog of classical particle flow 
along phase portrait lines. It reveals hidden features of quantum dynamics and extra complex- 
ity. Being constrained by conserved flow winding numbers, it also introduces topological order into 
quantum dynamics. 
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Phase portraits characterise the time evolution of dy- 
namical systems and are widely used in classical mechan- 
ics pp. For the conservative motion of a single particle, 
moving in one dimension x under the influence of a static 
smooth potential V(x) only, flow stagnation points deter- 
mine the topology of the overall classical Liouville flow 
in phase space [2] and are therefore of special interest. 
Situated on the x-axis wherever the potential is force- 
free (momentum p = and —dV/dx = 0), the local flow 
forms clockwise vortices around stagnation points at min- 
ima of the potential, maxima split the flow and therefore 
lie at the intersections of flow separatrices, saddle points 
of the potential lead to an elongated saddle flow pattern 
oriented along the x-axis. 

A single quantum particle, free to move in time t on a 
line parametrised by position x is described by a complex 
time-dependent amplitude function ijj(x;t). Wigner's 
function W(x,p;t) [3JH], with p the particle's momen- 
tum, is the quantum analog of classical phase space 
probability distributions. Structurally, W is a Fourier- 
transform of the off-diagonal coherences of the quantum 
system's density matrix g, i.e. 
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with Planck's constant h = h/(2ir). Unlike i/j or £>, the 
Wigner function only assumes real values, but these do 
become negative [3j |5], defying description in terms of 
classical probability theory |5 , thus revealing the quan- 
tum nature of a system [5H7J. 

For a single quantum particle, with mass m, moving in 
one dimension, the time evolution of W cast in the form 
of a flow field J(x,p;t) [8 , the 'Wigner flow' [9], has the 
two components 
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fulfilling Schrodinger's equation which takes the form of 
the continuity equation ^ + + ^ 



0. 
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Wigner flow has so far not been studied in detail [5HTT] . 

The dynamics of the harmonic potential, the most 
studied quantum case, e.g. of quantum optics [5j [7], 
amounts to a rigid rotation of the Wigner function 
around the origin of phase space. Only at the origin 
can a flow vortex form (Fig. 3 in Ref.[9]), just like in the 
classical case. The quantum harmonic oscillator and its 
isomorphism [12 , the free quantum particle [13], consti- 
tute exceptional, degenerate cases where lines of stagna- 
tion of Wigner flow occur, and not only isolated stag- 
nation points. This is due to the fact that for V oc x 2 
or V =const. we have in eq. Q J p = — just as 

in the 'classical limit' (h — >• 0) and consequently with 
W = we always find \J\ = 0. Because of these three 
facts (rigid rotation, classical form of J p and line forma- 
tion) the non-classical phase space features discussed here 
cannot be seen in the degenerate cases primarily studied 

so far Banna. 

Here we investigate the quantum dynamics of bound 
states of non-harmonic potentials, their classical coun- 
terparts exhibit regular flow [2], while their Wigner flow 
reveals rich non-classical features: 

Dependence of flow on the state of the system [14] 
leading to flow reversal [15. and directional deviation 
from classical trajectories [9], time-dependent quantum- 
displacement of classical stagnation points, occurrence of 
additional non-classical stagnation points and their dy- 
namics (even for conservative systems) , and conservation 
of the flow orientation winding number lj, see eq. be- 
low, carried by all Wigner flow stagnation points during 
all stages of their time evolution — including instances 
when they split from or merge with other stagnation 
points. 

The degeneracies of the degenerate cases leading to the 
formation of stagnation lines are lifted for non-harmonic 
potentials by the presence of terms with I > in J p and 
leads to the formation of separate stagnation points in- 
stead. The boundedness and continuity of wave functions 
of infinitely differentiable non-harmonic potentials and 
the unitarity of such systems' quantum dynamics induces 
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FIG. 1. (color online) Flow field around typical types of stag- 
nation points of Wigner flow with associated winding num- 
bers. This list is non-exhaustive. 



homotopies that keep all smooth changes of J in space 
and time around stagnation points smooth. We therefore 
expect that the topological structure of the Wigner flow 
field around stagnation points remains conserved. To 
monitor this we introduce the Wigner flow orientation 
winding number determined by the integral 
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along a closed (convex) loop C; here tp is the the orien- 
tation angle between the positive x-axis and the Wigner 
flow 'vectors' J. For 'empty' paths, not including a stag- 
nation point of J, uo = 0; for vortices lj = +1, see Fig [I] 
uj{C\t) is unchanged under continuous path deformations 
that do not push C across a stagnation point and as long 
as time evolution does not move a stagnation point across 
the loop. The winding number uj assumes integer val- 
ues only (assuming the integration path C does not run 
through a stagnation point) and is conserved. The sum 
of winding numbers of all stagnation points within a loop 
is conserved, even when they split or coalesce: the stag- 
nation points carry topological charge [IS] . 

To give an example, we concentrate on Caticha's [17] 
smooth, slightly asymmetric, double well potential 
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featuring high outer walls, and wells separated by a bar- 
rier of sufficient height, such that at least ground and first 
excited energy eigenstate tunnel through it, Fig.[2|a) in- 
set. To display all non-classical flow features listed above 
it suffices to investigate the balanced superposition 
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of ground ipo(x) = ipo cosh(x) exp [— (cosh 2 (x) + ax 
+ ^sinh(2x))] and first excited state ipi(x) = 
ipi [a + tanh(x)] ipo(x), with energies Eo and Eq + AE 
and normalisation constants and ipi , respectively [17] . 
The Wigner functions have to be determined numeri- 
cally [18 . Of the low energy states \I/ is the 'most dy- 
namic' in that it shifts all of the particle's probability 



distribution back and forth between left and right well. 
Its Wigner function W(t = T/4) at a quarter of the tun- 
nelling period time (T = 2irh/ AE) is displayed in Figs. [2] 
and 02 

The Wigner flow's state dependence has two aspects: 
the negativity of W reverses the flow, and the depen- 
dence of J p on W and V leads to sideways deviation of 
Wigner flow from classical phase portrait lines. In the 
case of eigenstates of the harmonic oscillator, the former 
leads to shear flow between neighbouring sectors of alter- 
nating polarity [9 . The latter deviation can remain mild 
for eigenstates of a weakly anharmonic potential [9], in 
our case it is very pronounced leading to the formation 
of non-classical vortices which are quantum-displaced off 
the x-axis and some of which spin anti-clockwise, Figs, [l] 
and 03 Indications of non-classical vortices seem to have 
been observed in 'chaotic' quantum systems before [T9] . 

Flow reversal affects quantum tunnelling. The wave 
function of a particle, tunnelling through a barrier, is 
coherent across it, implying that interference fringes of 
the Wigner function form in the tunnelling region; orien- 
tated parallel to the x-axis [6], Figs.^b) and [3) Neigh- 
bouring phase space regions contain strips of alternat- 
ing Wigner function polarity alternating their flow di- 
rection, Fig. [2^b) . The further apart the two wells, all 
else being equal, the finer the interference pattern [6] 




FIG. 2. (color online) Wigner function, Wigner flow and 
momentum distribution of state ^ (parameters H = l,m = 
1/2, a = 0.5, AE 0.5). (a) Wigner function W(x,p;T/4), 
projection to bottom shows its contours. Projections onto 
background walls show momentum and position probability 
distribution (blue and rose filled curves, respectively, in ar- 
bitrary units). Inset: Plot of Caticha-potential with wave 
functions for lowest two energy eigenstates shifted to their 
respective energy levels (dashed lines), (b) x-component of 
Wigner flow J x (Xs , p; t) at barrier top at position Xs ; bot- 
tom projection shows its contours (note its phase shift at 
t = T/2). Projections onto background walls show time and 
momentum projections (blue and rose filled curves, (J x (p))t = 
J T dr J x (X s ,p;r) and (J x (t)) P = f^dp J x (X s ,p;t), re- 
spectively, in arbitrary units). (J x (t)) p equals the probability 
current ^(^*d x ^ - ^d x ^*) at X s and features the sinu- 
soidal variation oc sin(27r£/T) expected of the tunnelling cur- 
rent of a 2-state system. 
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FIG. 3. (color online) Features of Wigner flow of an asym- 
metric double well potential (same parameters as in Fig. 1, 
at t = T/4). Colour-coded contour plot of Wigner func- 
tion W(x,p;T/4) (black contour lines) with an overlay of 
coloured arrows showing normalised Wigner flow J/\J\ (red 
arrows for W > and green for W < 0). The classical phase 
portrait is shown as a collection of thick white lines. All lo- 
cations of Wigner flow stagnation points are highlighted by 
symbols (cyan and red circles centre on clockwise and anti- 
clockwise vortices, respectively, yellow diamonds on separa- 
trix intersections, and the blue square on a p-directed saddle 
flow. The quantum-displacement of the vortices near the po- 
tential minima, towards the center, is clearly visible. 



and the stronger the resulting flow cancellation. Thus 
tunnelling can be described as a phase space transport 
phenomenon that plays out over large parts of quantum 
phase space [20] , frustrated through interference in phase 
space. 

For a low energy state, such as the positions of 
remnant vortices originating from the potential's minima 
positions are quantum-displaced inward, towards the po- 
tential barrier, Figs. [3] and |4j Different pictures emerge 
for other potentials [21] or systems in states more highly 
excited than Remnants of classical stagnation points 
can get substituted by more than one point each [21] and 
the overall flow pattern sensitively depends on the choice 
of state [2T] . 

The occurrence of stagnation points of Wigner flow 
away from the x-axis (p ^ 0) is impossible in classical 
physics. In the quantum case, even for conservative sys- 
tems, they can travel long distances in phase space and 
merge with or split from other stagnation points, Fig. [4] 

Their presence alters the topology of phase space flow. 

In the case of the Caticha potential in state we 
observe a string of vortices with alternating handedness 
aligned in the p-direction. They are located near the 
top of the tunnelling barrier, at position Xs = —0.258, 




FIG. 4. (color online) Wigner flow's stagnation points' posi- 
tions across phase space as a function of time. Same parame- 
ters as in Fig. [2] The tube surfaces show where in phase space 
the magnitude of Wigner flow is small J(x,p;t) 2 = 3 • 10 -5 . 
They are displayed over 120% of one oscillation period (t — 
—0.1 • T. . . 1.1 • T); the rainbow spectrum is matched to T, 
red-orange for t — 0, via green, cyan at t — 0.5T, through 
blue and purple back to red-orange. Because of the period- 
icity of the two- state scenario, the red-orange-yellow torus is 
seen twice at beginning and end of the time window overlap. 
At the core of all tubes lie time-lines of stagnation points. 
Their movement through phase space leads to their mergers 
and splits. The grey flow integration lines are guiding the eye 
past vortices and separatrices, they do not represent physical 
flow since they integrate J(x, p; t = 0.05) at fixed time. The in- 
coordinate is shown from left well minimum at Xl — —2.095 
to right well minimum at Xr — 1.514, the position of the 
associated vortices (©in left well and ® in right), at p = 
and just inside the plot region, confirms the inward quantum- 
displacement of these stagnation points. The remnant of the 
classical separatrix stagnation point (p = 0, labelled ©) at 
position Xs + 6xs(t) ~ —0.3 follows a slightly bent path. 
It forms the handle of 'Poseidon's green trident', which turns 
blue as time progresses and then the remnant separatrix stag- 
nation point merges with an incoming counterclockwise vor- 
tex in the purple wrist of the four-fingered claw, ceasing to be 
a separatrix intersection. Shortly after, a separatrix intersec- 
tion temporarily emerges in conjunction with a vortex to form 
the torus, subsequently we are 'back to Poseidon's trident'. 
Mergers or splits are constrained by topological charge [TB] 
conservation. 



pinned to the lowest points (in the p-direction) of the 
Wigner function's interference pattern's meniscus-shaped 
zero-lines, Fig. [3j Over time they travel downward with 
the downward movement (in the p-direction) of the in- 
terference fringes, Fig. [4] When reaching the x-axis they 
coalesce with the remnant separatrix intersection point. 
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Upon their merger the local flow switches to a saddle 
point's flow pointing in the p-direction, which is impos- 
sible in classical phase space, Figs. [I] and [3] 

Compared to classical phase portraits, Wigner flow 
reveals non-classical features and added complexity of 
quantum phase space dynamics; at the same time it pro- 
vides, through the conservation of cj, the basis for an 
analysis of its topological invariants, ordering this com- 
plexity. Our numerical findings [18] confirm that winding 
number conservation applies to isolated stagnation points 
as well as those that interact with others through splits 
or mergers, Fig. [4] As an example one can create an inte- 
gration loop C around the torus of Fig. |4j since the torus 



is formed of a vortex and a separatrix crossing its overall 
charge is zero, as required given that it appeared out of 
'nowhere'. 

We conjecture that all stagnation points of the Wigner 
flow of all superpositions of bound states of non-harmonic 
potentials carry conserved topological charge [TS] . 

Systems that have been studied using quantum phase 
space techniques can be analysed using Wigner flow. 
Such systems arise for example in chemical quantum 
dynamics [8] [TO] [22], 'non-linear' quantum processes in 
closed single particle [11 or open multi-particle [23H25] 
systems, classical electromagnetic fields [26] and multi- 
band semiconductor physics [27] . 
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